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Abstract. We report theoretical values for the unscreened plasma frequencies Qp of 
several Fe pnictides obtained from density functional theory (DFT) based calculations 
within the local density approximation (LDA) and compare them with experimental 
plasma frequencies obtained from reflectivity measurements on both polycrystalline 
samples and single crystals. The sizable renormalization observed for all considered 
compounds points to the presence of significant many-body effects beyond the LDA. 
From the large values of the empirical background dielectric constant Soo ~ IS- 
IS derived from reflectivity data, we estimate a large arsenic polarizability a^s ~ 
9.5 ±1.2 A'^ where the details depend on the polarizabilities of the remaining ions 
taken from the literature. This large polarizability can significantly reduce the value 
of the Coulomb repulsion Ud ^^4 eV on iron known from iron oxides to a level of 2 eV 
or below. In general, independently on such details, this result points to rather strong 
polaronic effects as suggested by G.A. Sawatzky et ai, in references arA"w:0808.1390 
[TT] and arA"w:0811.0214p^. Possible consequences for the conditions of a formation 
of bipolarons are discussed, too. From the extrapolated fiSK (muon spin rotation) 
penetration depth data at very low-temperature and the experimental value of the 
unscreened plasma frequency we estimate the total coupling constant Xtot for the 
electron-boson interaction within the framework of the Eliashberg-theory adopting 
an effective single band approximation. For LaFeAsOo.gFg.i a weak to intermediately 
strong coupling regime and a quasi-clean limit behaviour arc found. For a pronounced 
multiband case we obtain a constraint for various intraband coupling constants which 
in principle allows for a sizable strong coupling in bands with either slow electrons or 
holes. 
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1. Introduction 

Naturally, shortly after the discovery [1] of superconductivity at relatively high transition 
temperatures up to 56 K in several iron pnictide compounds despite all efforts world 
wide, the underlying pairing mechanism and many basic physical properties both in 
the superconducting and in the normal state are still debated. This concerns also the 
nonsuper conducting parent compounds, which often show magnetic phase transitions. 
In such a situation a combined theoretical and experimental study of selected systems 
as well as a detailed comparison of various members of the fastly growing family among 
themselves and also with other more or less exotic superconductors, especially with the 
high temperature superconducting cuprates (HTSC) might provide useful insight into 
related problems. In this context the unscreened plasma frequency Qp determined from 
optical measurements or electron energy-loss spectroscopy (EELS) at frequencies near 
the visible range 

= 1^ (1) 

m* 

is a physical quantity of central interest in the present paper for several reasons, (i) It 
provides direct insight in the dynamics of free charge carriers that are responsible for 
the superconductivity. In particular, it reflects also the possible influence of many-body 
effects in changing their effective mass m* and measures the density of paired electrons 
(condensate or superfluid density), if as in the BCS-case n = Ug holds, (ii) According 



New insight into the physics of iron 



3 



to the famous Uemura plot [2] for T^. vs. Ug/m*, simple linear relations between the 
two quantities may be expected for many exotic type-II superconductors, and Vt^i would 
then be a direct measure of Tc. Such plots have been presented both for underdoped 
HTSC and for the pnictide superconductors as well [2l[3llll[5l[6l[7l[8]. In some cases, 
especially for members of the 122-family, i.e. for Ba(Sr)Fe2As2 derived superconductors 
its applicability has been put in question [8|, [9]. Anyhow, since a Bose condensation 
of bipolarons is naturally directly related to the superfluid density [10], bipolaronic 
scenarios such as the recently proposed one by G.A. Sawatzky et al. [HI [12] are of 
special interest. 

The unscreened plasma frequency is also of considerable interest for the analysis of 
resistivity (p(T)) data analyzed within the frequently used Drude model. In particular, 
the high-temperature linear T-region is often used to extract the transport electron- 
phonon (el-ph) coupling constant Xtr [13] using the simple relation 

A - (2) 

For standard superconductors Xtr does not differ much from the Fermi surface averaged 
el-ph coupling constant Xei-ph which governs the transition temperature Tc of the 
superconducting state and it is also responsible for the renormalization of the electronic 
specific heat of a Fermi liquid at low-temperature. Using our previous very first estimate 
of Op for LaOo.gFo.iFeAs [H], equation (2) has been used [131 [15] to consider a classical 
el-ph scenario with A ~1.3 for the high Tg- values and for a high-temperature saturation- 
like behaviour of the resistivity in the Pr based 1111-iron pnictides. However, the 
present improved knowledge and data do not provide support for such a point of view 
(see below). For more realistic estimates within a two-band model approach see reference 

Finally, Qp and A determine also the pair-breaking parameter /? = QpPQ/8n{l + X)Tco 
for various unconventional pairing states [T71 jlHl [19], |20] : 

-u^).jum)-^'a). (3) 



where ip{x) is the digamma function and T^o denotes the transition temperature being 
at maximum in the clean-limit where the residual resistivity po ~^ 0- 

Here we report on the optical analysis of superconducting and nonsuperconducting 
parent compounds. Besides our own data also published data of other groups will 
be discussed and used in our analysis. We briefly introduce the employed codes for 
the determination of the electronic structure and the unscreened plasma frequencies, 
which we regard as quantities not affected by several possible many-body effects such 
as correlation, nonadiabatic, polaronic and polarization effects being under debate at 
present. With the aid of penetration depth data near T = taken from pSR (muon spin 
rotation (relaxation) )-data we are in a position to estimate the total coupling strength 
between the charge carriers (both electrons and holes) and some still unspecified bosonic 
mode(s). The analysis of the background dielectric constant Eoo in the range of the 
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screened plasma frequency as derived from optical measurements sheds light onto the 
role played by correlation and polarization effects which is a central problem of general 
interest in modern solid state physics independently on the peculiar physics of high- 
temperature superconductivity. 

2. Reflectivity and its analysis 

2.1. Studies on polycrystalline LaOi-^^F^-FeAs samples 

Polycrystalline samples of LaOo.gFo.iFeAs [211 122] as well as of the undoped parent 
system LaOFeAs (also called (1111) compound) have been prepared as described in 
previous work [5], [23]. In general, the situation here parallels that for the HTSC. Right 
after their discovery, also only polycrystalline samples were available, and have been 
studied by reflectance measurements. Further, also the HTSC are quasi-2D metals with 
the highest plasma energy (PE) value for light polarization within the (a,b) crystal plane. 
In the case of HTSC, there is a number of publications devoted to polycrystalline pellets, 
exactly as in our studies of La(0,F)FeAs, which demonstrate that the reflectivity onset 
seen as a kink in the measured curves is a good measure of the in-plane PE value (see 
e.g. Kamaras K et al. [21]). Such studies already provided very important information 
on the optical properties of HTSC, e.g. they were used to study the doping dependence 
of the in-plane PE, and most importantly these results later on were corroborated by 
single crystal studies (see e.g. reference [25]). 

For the optical measurements reported here the pellets have been polished to 
obtain appropriate surfaces. The reflectance measurements have been performed using 
a combination of Bruker IFS113v/IFS88 spectrometers. This allows us to determine 
the reflectivity from the far infrared up to the visible spectral region in an energy range 
from 0.009 up to 3 eV. The measurements have been carried out with different spectral 
resolutions depending on the energy range, these vary form 0.06 to 2 meV. Since the 
observed spectral features are significantly broader than these values, this slight variation 
in resolution does not impact our data analysis. All reflectance measurements were 
performed at T = 300 K. 

In figure 1 we show the reflectance of our LaOi_2;F^FeAs samples {x = 0,0.1). 
The reflectance significantly drops between 0.13 and 0.4 eV, thereafter it shows a small 
increase until about 0.6 eV before it starts to slightly decrease again. In addition, at 
about 2 eV another small upturn is visible for LaOo.gFo.iFeAs. We attribute the two 
small features at about 0.6 and 2 eV to weak electronic interband transition at the 
corresponding energies. This observation is in accord with the behavior of the optical 
conductivity cr(co') predicted by Haule et al. [26] using dynamical mean-field theory 
(DMFT), where As 4p to Fe 3d interband transition slightly below 2 eV and a weak 
feature near 0.6 eV for the undoped parent compound LaOFeAs have been found. Again, 
our measurements on powder samples do not allow the extraction of their polarization. 
The steepest edge centered at about 250 meV with an high energy onset at about 
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Figure 1. Reflectance of polished undoped LaOFeAs, optimally doped 
LaOo.gFo.iFeAs, as well as an arsenic-deficient sample LaOg.gFo.iFeAsi-^ with (5 sa 0.1. 
All three curves are characterized by a rather pronounced kink which is attributed to 
a plasma edge at about 380 meV and 395 meV, respectively. For the interpretation 
of the broad edge-like feature near 240 meV in the special case of LaOo.gFg.iFeAsi-i 
see figure 2 and text. In addition, there are signatures of interband excitations around 
0.6 eV (and for LaOo.gFo.iFeAs at about 2 eV). 

380(20) meV for LaOFeAs and 395(20) meV for LaOo.gFo.iFeAs represents the plasma 
edge or plasma energy, which is observed at relatively low values. In consideration of the 
quasi-2D character of the electronic states in LaOi^^F^FeAs and the expected strong 
anisotropy of the plasma energy for light polarized || to the (a,b) and c crystal axes, 
respectively, the onset of the plasma edge in figure 1 gives a reasonable value of the 
in-plane value (i.e. for light polarized within the (a,b) crystal plane). Consequently, the 
in-plane plasma energies for the two LaOi_a;Fa,FeAs samples (undoped and 10% electron 
doped) is almost independent of the doping levels studied here. 

Within a simple Drude model, the value for the plasma energy allows a first, rather 
crude estimate of the charge carrier density in LaOi_a;Fa;FeAs. Within this model, 
the observed screened plasma energy, cUp, is given by oj^ = Anne^/meoo, with n, e. 
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Figure 2. Idealized in-plane Drude-Lorentz model for the loss function (left) and 
the reflectivity (right) of optimally doped LaOo.gFo.iFeAs, as well as an arsenic- 
deficient sample LaOo.gFo.iFeAsi-^ with S « 0.1. For a qualitative comparison with 
experimental polycrystalline data see figure 1. For fit parameters see text. 



m, and Eoo representing the charge carrier density, the elementary charge, the effective 
mass of the charge carriers and the background dielectric screening due to higher lying 
electronic excitations [27] . Adopting an effective mass equal to the free electron mass, 
and ^oo ~ 12 [2H], one would arrive at a rather low charge density of about 9-10"^'' cm~^. 
We note, that this estimate ignores the strongly anisotropic nature of the electronic 
orbitals near the Fermi energy in LaOi_^F^FeAs. A more detailed discussion of the 
unscreened in-plane plasma energies found well above 1.3 eV but below 2 eV in all 
pnictide superconductors is given in the next sections see also table 2. 

Finally, we examined also an arsenic-deficient sample LaOo.gFo.iFeAsi.^ with 
6 ~ 0.1 [22l [29]. Here at first glance a plasma-edge like feature near Up = 0.24 eV, 
only, would suggest that the screened plasma frequency is considerably lowered, (see 
figure 1). Adopting a linear scaling of the background polarizability ^oo ~ 1 with the 
As-concentration, one arrives at ^oo ~ 10.5 instead of 12 for the non-deficient sample. 
Within a nonlinear scaling suggested by the Clausius-Mossotti relationship between the 
dielectric constant and the ionic polarizibilies discussed in more detail in section 4.2 one 
estimates even a slightly smaller value of ^oo ~ 9.97. Thus, one would arrive at Qp = 0.79 
to 0.76 eV for the unscreened plasma frequency. However, such a strong change of Qp 
by about 0.6 eV (!) would imply very large changes in the electronic structure (e.g. by 
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filling the hole pockets and an additional mass enhancement) or alternatively a twice 
as large ^oo- Both scenarios are difficult to imagine microscopically. For this reason we 
will consider a third possible explanation with an almost unchanged VLp and a slightly 
reduced Soo ~ 10. We ascribe the clearly changed optical properties to the presence of 
arsenic derived electrons and adopt that an absorption peak described by a Lorentzian 
curve like that absorption near 0.66 eV in the BaFe2As2 single crystal data with about 
5% arsenic vacancies [30] (see below) might be mainly responsible for that unexpected 
feature. Simulating qualitatively the polycrystalline reflectance data shown in figure 
1 by an effective Drude-Lorentz model we arrive at a corresponding Lorentzian curve 
peaked near 1.16 eV and a slightly lowered unscreened plasma frequency Vtp =1.3 eV 
(A quantitative fit would require a more or less realistic account of the shapes of the 
microcrystalline grains and their anisotropic dielectric properties which, however, is 
without the scope of the present paper.). The oscillator strength of this Lorentzian- 
peak is expected to be about twice as large which is supported by our fit value of 6.5 
eV, see figures 2-4 and compare table 1. The different binding energy might be related 
to the rather different As-environment in the 1111 and the 122 compounds with no 
La(0,F) layer in the latter and an enhanced pinning provided by the La'^"'" ions being 
close to the As-sites whereas in the former case it is provided by the less effective 
pinning due to Ba^'''-ions. The resulting loss function and the optical conductivity 
are shown in figure HI The obtained optical conductivity resembles the mid-infrared 
absorption observed by Yang et al. in a Bao.55Ko.45Fe2Asi_5 single crystal [31], [30] with 
a Tc ~ 30 K slightly above our arsenic-deficient La-Ill sample. This observation gives 
further support for our vacancy derived bound-electron scenario. In this context it 
would be interesting to study analogous effects in the related superconducting FeSeo.ss 
and other substoichiometric iron telluride and selenide compounds |32j . 

Here we would like to mention that a very small value for the unscreened plasma 
frequency has been derived from an ellipsometry study of LaOo.gFo.iFeAs and a 
subsequent analysis applying the effective medium approximation (EMA) to obtain 
a^s{u) for a polycrystalline sample [2H|. For spherical grains and aab ^ cTc as expected 
for a quasi-2D metal (as is the case for LaOo.gFo.iFeAs), the EMA predicts (7^^ ~ O.ScXaft, 
which then would correspond to a slightly larger value of Qf^ = 0.86 eV. In addition 
some spectral weight from the Drude part was lost in reference [28] due to the fitting 
of broad Lorentzian-peaks also for the two low-a; interband transitions ai and a2 at 
0.75 and 1.15 eV, respectively. Adding this lost weight to the Drude term we arrive 
at an enhanced corrected Drude weight by a factor of 1.28 which corresponds to a 
corrected fi"* ~0.975 eV. Moreover, since this value has been derived mainly from the 
low-u, damped Drude region (below 25 meV), it is still renormalized by electron-boson 
coupling, and a high-energy boson well above 25 meV has been assumed to explain the 
high Tc at weak coupling. Using a weak total electron- bo son coupling constant A ~ 0.61 
( see section 5.1) one estimates Qp =1.24 eV, close to a value of 1.37 eV suggested from 
the plasma energies in figure 1 and dielectric screening (see above). 
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2.2. BaFe2As2 single crystal studies 

With the successful synthesis of single crystals, more detailed investigations of the 
optical response of iron pnictide superconductors become possible. In particular, 
the background dielectric constant ^oo can be determined more safely. Reflectivity 
data of Bai_a;Ka;Fe2As2 and SrFe2As2 single crystals have been reported recently, and 
their analysis corroborates what has been derived from the polycrystalline samples as 
discussed above. In figure 3 we present reflectivity data of a BaFe2As2 single crystal 
and a light polarization within the (a,b) crystal plane. The investigated single crystal 
with the dimension of about 0.5x0.5x0.1 mm was grown from a tin flux. Elemental Ba, 
Fe, As were added to Sn in the ratio 1 : 2 : 2 : 15 and placed in a corundum crucible 
which was sealed in an argon filled silica ampoule. The sealed ampoule was heated in 
steps to 850° C with a rate of 20 K/h, hold for 36 hours and cooled down. The tin flux 
was removed by solving in diluted HCl. The resulting gray, crystals with metallic luster 
were washed in distilled water several times and finally dried. 

The studied single crystal shows a thin platelet-like shape and a shiny (0, 0, 
1) surface. The thorough characterizations via specific heat, resistivity, magnetic 
susceptibility and chemical composition analysis have been carried out on the crystals 
obtained from the same batch. Similar physical properties to those in the crystals 
reported in references [30] have been obtained. The EDX composition analysis of several 
specimens of this batch resulted an average composition of Bao.95Sno,o5Fe2As2 and we 
expect ~ 5% Sn to be incorporated in the crystal structure. For the consequences of 
an As-deficiency with respect to low-frequency optical properties see below and for the 
upper critical field Bc2{T) of the K-doped system see [29] . 

Also these data are characterized by a "kink" -like feature (but not by a sharp 
minimum as expected in the case of a clean single crystal!) in the reflectivity around 
0.4 eV, with a steeper decrease at lower energy values. The polycrystalline samples as 
analyzed in the previous section did not allow a quantitative analysis of the reflectivity 
values due to relatively rough surfaces, but in the case of single crystal data as shown 
in figure 2 this is feasible, and figure 2 also shows a fit of our reflectivity data within a 
Drude-Lorentz model: 

n2 3 ^2 

eiuj) = e^- ^ +E- 1^^- (4) 

This fit describes the data very well and allows the extraction of more physical 
information. The resulting fit parameters are summarized in table 1. 

One important finding of our fit is the value for the screened plasma frequency 
Up = Qp/ ^/e^ of about 400 meV, which is in very good correspondence to what has 
been derived from the polycrystalline LaOi_a;FxFeAs samples above. This demonstrates 
that this value is rather independent of the stoichiometry but is a characteristic feature 
of the FeAs planes common to all FeAs based superconductors. Furthermore, in order 
to obtain a reasonable description of our data within the Drude-Lorentz approach, it is 
necessary to include Lorentz oscillators above and below this unscreened plasma energy 
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Energy (eV) 

Figure 3. Reflectivity of a BaFe2As2 single crystal and a light polarization within 
the (a,b) crystal plane (open circles). The solid line represents the result of our Drude 
Lorentz analysis (see text). 

(see table 1) at least for the BaFe2As2 (also called (122) compound), and other (122) 
systems. 




Figure 4. Experimental in-planc loss function (left) and optical conductivity (right) 
of a BaFc2As2 single crystal and the same for the suggested minimum Drude-Lorentz 
model for the arsenic-deficient La-1111 system. Adopting wq,! =1.16 cV 70,1 = 3.3 eV 
and rip4=6.5 cV for the "bound" electrons as well as Fp =0.03 eV and fip =1.3 eV for 
the free electrons for the latter. Compare table [T] and equation 
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Table 1. Parameters derived from a Drude-Lorentz fit of the reflectivity of a BaFe2As2 
single crystal. Here the following notation has been used: ^Ip = Drude plasma 
frequency. F = Drude width, wo.j = Lorentz oscillator position, 7^ ~ Lorentz oscillator 
width, ujpj = individual Lorentz plasma frequency ( or oscillator strength), and eoo = 
background dielectric constant). 



j 




7, (eV) 


ujpjieV) 


r, (eV) 




1 


0.2 


0.16 


0.79 


0.031 


1.58 14.9 


2 


0.66 


2.9 


4.31 


0.087 [33] 


1.6 [33] 


3 


0.93 


1.63 


1.86 







Regarding the clean (1111) systems, so far clearly visible interband transitions 
have been observed only above the plasma edge [211 128] . An inspection of the atomic 
occupancies given for such single crystals obtained from the Sn-fiux by [30], reveals 
a non-neglible amount of disorder, first of all of arsenic-vacancies about half as large 
as in the polycrystalline arsenic-deficient sample reported above. For this reason we 
assume that the low-frequency Lorentzian-peaks describe in both cases transitions from 
As-vacancy introduced, more or less localized, states below the Fermi-level Ep but well 
above the usual As 4p states which occur about 2 eV below Ep. These states are 
thought to accomodate the large amount of electrons compared for instance with its 
free electron Hall number uh ^ 5 ■ 10^°cm~^ [29]: ~ CyacZ /Vceu ~ lO^^cm"'^ which would 
otherwise heavily overdope the La-1111 system or electron dope the nominal 122 parent 
compound, where Z <3 is the effective anionic As-charge, c^ac ~ 0.05 to 0.1 denotes the 
atomic As-vacancy concentration an v^^n =70.75A^ denotes the molar cell volume. 

Although the reflectivity curves for FeAs based single crystals as presented in 
figure 3 and in the literature so far [31], [331 EH ES, [361 [31] do not agree in all aspects, 
there are essential parameters in common, which give quite a lot of insight into the 
physics of these systems. First, interband excitations are present in the energy ranges 
around 0.2-0.3 eV and 0.6 - 2 eV which most likely arise from the manifold of iron related 
bands that cross the Fermi level [HH]- Second, the background dielectric screening in 
the energy range of the plasma energy is quite large with values of about 12 to 16. 
Noteworthy, comparable values have been already observed occasionally also for other 
transition metal compounds [39] . Third, the value of the (unscreened) plasma frequency 
is in the range of 1.4 - 2 eV, rather independent of the doping level, and importantly, 
it is substantially smaller than what is predicted from density functional (DFT) based 
band structure calculations. The latter issue will be discussed in detail in the following 
sections. 

3. FPLO analysis of the unscreened plasma frequency 

The theoretical plasma frequencies have been taken from DFT based calculations in the 
local (spin) density approximation (L(S)DA). We used the full potential local orbital 
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code [10] in version 7 with the standard double numerical basis. The theoretical 
unscreened (Drude) plasma frequencies are derived from the Lindhard intra band 
expression Vtp (see references [SI US] for the isotropic case and equation ^ below 
for the straightforward anisotropic generalization). They basically represent the Fermi 
surface average of the Fermi velocity squared. 

It has been checked that the results are converged within a few percent with 
respect to the /c-integration scheme. We used the experimental lattice constants and the 
theoretical arsenic position. However, the frequencies seem not to be too sensitive to the 
lattice parameters. In the case of doped systems the simple virtual crystal approximation 
has been used. Thus, here we report calculations of the eigenvalues of the anisotropic 
conductivity (dielectric function) tensor 

^''^li,.o. = ^^e''N{Ep)vl a = x,y,z (5) 

where the Fermi velocity components Va have been averaged over all Fermi surface sheets 
(FSS) and NiEp) denotes the total density of states at the Fermi level Ep for both spin 
directions. In the tetragonal phase the 'xx' and 'y?/' component do coincide and we 
will them denote as in-plane plasma 'a6'-frequency whereas the out-of plane 'c'-plasma 
frequency is given by the ''zz' component. In terms of individual Fermi surface sheets 
(bands) equation ([5]) can be rewritten as 

K,o. = Y.K, (6) 

i 

where i denotes the band index. In the pnictides under consideration all DFT- 
calculations predict up to five bands (FSS) at Ep consisting of electron (e/) and hole 
ih) FSS. For the sake of simplicity the discussion of multiband effects all these FSS can 
be lumped into one effective e/-FSS and one /i-FSS, thus resulting in frequently used 
two-band models. 

Due to the layered structure of the pnictides one has fi^^ ^ Jl^. The larger in- 
plane distance and the smaller transfer integrals in c-direction for the (1111) systems as 
compared for instance with the (122) systems is mainly responsible for the relative small 
values of both in-plane and out-of plane comnponents. For comparison we included also 
calculations by Singh [13] for the case LaOFeP and Boeri et al. [H] for the case of 
LaOFeAs and LaOFeP where different bandstructure codes have been used. We note 
that different codes results in essentially the same unscreened plasma frequency. All 
are clearly a factor of 1.5 to 2 larger compared with the experimental values (see table 
2), pointing to the presence of many-body effects which for instance may cause a mass 
renormalization. Several possibilities to resolve this issue will be discussed below. 

The inspection of the calculated and empirically found unscreened plasma 
frequencies shows a very sensitive dependence on the Fe-Fe- distance which in the 
tetragonal phase is given by the in-plane lattice constant a/\/2. Since fip oc oc td^a, 
where t^^ is a representative value for the in-plane Fe-Fe 3(i transfer integrals, a simple 
power-law has been predicted considering the related Slater-Koster integrals [15]. As a 
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Table 2. Comparison of calculated (DFT) and experimental unscreened in-plane 
(a6) plasma frequencies measured in eV for polycrystalline samples (PC) and single 
crystals (SC). For completeness, the theoretical out-of plane values are shown, too. 
In the references the first citations notation denotes the experimental optical data (if 
available) and the second one the structural data employed in the DFT calculations. 
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1.72 
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1.58 


SC 


present work 








1.39 


SC 


[36] 


K^Bai_a;Fe2As2 


2.63 


0.8 


1.5 


SC 




Ko.45Bao.55Fe2As2 


2.63 


1.72 


1.9 ± 0.2 


SC 


[31] 


CaFe2As2 


2.95 


1.96 








LiFeAs 


2.9 


0.83 






m 


NaFeAs 


2.68 


1.16 






m 


LaOFeP 


2.37 


0.47 


1.85 


SC 


m 


result we expect 


^p,o/(a/ao) 


U 

1 






(7) 



where ao =4 A has been introduced for convenience. According to Harrison [15] we 
adopt 1/ = 4 and arrive at at f2p^o = 1-4 eV. Applying the approximate expression ([7]) to 
rare earth-1111 systems, where to the best of our knowledge no high-frequency optical 
data are available, and considering for instance two superconductors with the highest 
transition temperatures Tc > 50 K: SmOo.ssFeAs and NdOo.ssFeAs with a = 3.897 A 
and a = 3.943 A we estimate flp = 1.55 eV and 1.48 eV, respectively. 

Having now a large enough set of available empirical plasma frequencies f2p ~ 1.4 
to 1.7 eV we are in a position to re-estimate the mentioned above transport coupling 
constants [13], |15]. With these values using equation (2) we would arrive, instead at 
X Xtr = 1.3, at a superstrong coupling regime A ~ 4 to 5 (!). Then we should have to 
ask, why is Tc of the iron pnictides so low? In our opinion other scattering mechanisms 
than the usual e/-j>/i interaction, in the el-eloi the el-spin fluctuation channels, contribute 
to p{T) and possibly dominate it in the quasilinear region. Adopting a weak el-ph 



New insight into the physics of iron 



13 



coupling constant ~0.2 like the one estimated by Boeri et al. [H] a second linear region 
might be achieved beyond 600 K, only. To the best of our knowledge there are no 
experimental data available at present for that region. 

4. Discussion of possible related many-body effects 

4-1. The Mott-Hubbard scenario 

Our LDA-FPLO band structure calculations for the nonmagnetic ground state of 
LaOFeAs provide a value for the plasma frequency within the (a,b) crystal plane of 
flp^^ = 2.1 eV and a small value for a polarization along the c axis of 0.34 eV in 
accord with 2.3 eV and 0.32 eV given in reference [B]. Note that these "bare" values 
do not describe the screening through e^o caused by interband transitions. In order 
to get the measured values of Up = 0.395 eV for the fluorine doped and 0.38 eV for 
the undoped parent system an unusually large value of Eoo ~ 28.3 to 33.9 would be 
required which seems to be unrealistic and instead intermediate values Eoo ~ 12 to 15 
considerably larger than those for cuprates (4 to 6) would be expected [28]. Then, for 
the empirical unscreened plasma energy Qp a value about 1.4 eV would be expected 
on the basis of our reflectance data. Thus the question arises, what is the origin of 
these renormalizations? In view of the transition metal iron it is natural to consider 
moderate short range Coulomb interaction induced correlations as a possible reason. 
To get more insight into the role of possible correlation effects adopting the parameter 
set of Haule et al. [26] we calculated at first numerically the interband contribution to 
the complex dielectric function e(a;) from the in-plane DMFT optical conductivity cr(a;) 
between 1 and 6 eV for the undoped LaOFeAs given in figure 5 of reference [26]. From 
its static value we obtained Eoo = 5.4, only (see figure 5). Taking into account the small 
differences in the optical spectra between the undoped LaOFeAs and the optimally 
doped LaOo.gFo.iFeAs as shown in sections 2 and 3, we will not distinguish between 
their ^oo-values in the following discussion. (Using the ionic polarizabilities of F"^ and 
0~^ given in section 4.2 and the unit cell volume of f =142.32 [1] for LaOFeAs, one 
estimates from equation (9) a slight change of its Eoo to 11.61 or 11.25 for the parameter 
sets of references [HI |19] and [50], [511 [521 [531 [SH [55], respectively, compared with 11.5 
to 12 for LaOo.gFo.iFeAs [28].) In the strongly correlated (narrow band limit) limit 
Ud ^ W, for a very qualitative discussion, Eoo can be written approximately in the 
following way 



where W ~2 eV denotes the total Fe 3d bandwidth, Epd ~ 2eV is the oscillator strength 
for the arsenic-iron Ap-FeSd interband transitions and fdd describes the transitions 
through the Mott-gap. From our numerical simulations for Ud =4 eV we estimate 
fpd = 4.6 eV and fdd =8.4 eV. The resulting curve is shown in figure [S] (left). 

Despite the crude estimate of equation [8] in the weak [/-limit, in general for smaller 
values for the onsite Coulomb repulsion Ud within the DMFT-approach would result in 



oo 



1 + 




(8) 
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Figure 5. Intcrband contribution to the real part of tlic dielectrie funetion £i{uj) of 
LaOFeAs using the DMFT-calculation for its iy{Lo) as given by Haule et al. in reference 
in figure 5 therein with Ud = 4 eV. 



somewhat larger values due to an energy shift of the interband transition from the 
lower Hubbard band. Hence, a corresponding systematic DMFT study with sufficiently 
smaller Ud and J^-values would be of interest. Alternatively to figure 6 (left) one might 
expect that for small L'^-values there should be no visible infiuence of U at all. Then a 
negative curvature instead of the positive curvature suggested by equation (8) might be 
realized. A corresponding guess taken from our previous paper [21] is shown in figure 6 
(right). Anyhow, a more sophisticated study should shed light on the "screening" 
mechanism of the bare Coulomb repulsion. In this context it is noteworthy that the 
arsenic-polarization (solvation) mechanism proposed by Sawatzky et al. [Til [12] briefly 
discussed in the next subsection 5.1.2 seems to give a similar behaviour as equation (8) 
(see figure (7) (right)). 

In contrast to Haule et al. [26], Anisimov and Shorikov et al. [561 [57] argue 
that LaOFeAs is in an intermediate U regime but strongly affected by the value of 
the Hund's rule exchange Jd- Hence, ^oo should be strongly increase and the adopted 
relatively strong Coulomb repulsion [/^ = 4 eV is not compatible with the available 
optical data implying a much larger > 10-15. Thus, we are obliged to consider a 
weakly correlated scenario. 

A somewhat similar attempt to describe the experimental reflectivity R{uj), a 
reduced plasma frequency and "bad metal" effects, i.e. strong damping effects (flrst 
considered by Haule [26] like in a marginal Fermi liquid ) has been undertaken in a 
recent study by Laad et al. [58] using a simplifled version of the dynamical mean- 
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field theory. Although at first glance a seemingly similar behaviour has been obtained, 
quantitatively the results are not very convincing because (i) using = A eV and at 
the same time ignoring the hybridization with As 4p states reduces the minimal eight- 
band model to an effective five-band model where reduced effective Coulomb interactions 
and in particular also the onsite term [/-values should be considerably reduced, (ii) the 
important screening by the background due to pd and other non-dd interband transition 
described implicitly also by in our approach has been ignored in calculating the 
reflectivity R{uj), (iii) a comparison with results from measurements on powder samples 
requires a proper treatment of the anisotropic grains, for instance, at least approximately 
within an effective medium theory. 

4.2. The weakly correlated As-polarization scenario 

Taking into account the criticism of strong on-site Coulomb repulsion Ua and Hund's rule 
coupling exchange Jd mentioned above, the recently proposed novel, quite interesting 
polarization (solvation) scenario [11] for the pnictide components As, Te, Se, or P as 
seen by the Fe 3d electrons is of special interest. In this context we will compare our 
empirical values of Eoo with the contribution of the polarizability from As alone, a~As, 
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= 4 eV (DMFT) 
Experiment 
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Figure 6. Semi-empirical relation between the total coupling constant from the mass 
enhancement entering the penetration depth vs. interband screening modeled by the 
dielectric background constant using equation (8) (left) and a first estimate proposed 
in reference [21] (right) (see also figure 7 (right)). 
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as adopted by Sawatzky et al. [H]. For this purpose to get a first rough estimate, we 
employ the well-known Clausius-Mossotti (also known as Lorenz-Lorentz relation) for 
the sake of simplicity ignoring anisotropic screening effects 

atot = T ^> " = ^ = 1 1 TTT - 2 > 1, (9) 

+ 2 3 l-batot/iv) 

(Note that the Lorentz factor 6 = 47r/3 has been derived for the case of cubic symmetry.) 
Using for LaOo.gFo.iFeAs £oo=12 and v = 140. SA^ with two As-sites per unit cell we 
arrive at atot = 26.35 and atot = 13.17 A'^ per As-site already close to 9 to 12 A^ 
adopted by Sawatzky et al. [H]. For BaFe2As2 with £:oo=14.9, v = 203. 96A^ with four 
As sites per unit cell one obtains atot =40.05 A^ which corresponds to atot = 10.01 A^ 
per As-site. Since the other atomic sites are expected also to contribute to atot, we briefly 
consider their influence adopting the simple expression provided by the polarizability 
additivity rule [H] 




aAs = atot - TT yi"ion, /(I - Cy^As), (10) 



where the sum runs over all remaining ions(atoms), Nc is the number of As-sites per unit 
cell, and Cy^As ~ 0.1 denotes the atomic concentration of As- vacancies. Then it is clear 
that these ignored ionic contributions may somewhat change our empirical estimates 
of the dominant As-contribution we are mainly interested in. In fact, adopting for the 
contribution of Ba'^+ 1.7 A^ and -1.0 A^ for that of Fe2+ [50l [55] one arrives at a slightly 
larger value for ckas ~ 10.32 A^, and a smaller value of ckas^- ~ 7.89 A^ is obtained using 
instead the parameters ap^2+ = 1.339 A^ and 03^2+ = 1-57 as recommended in reference 
[18] . Taking into account the presence of As- vacancies in the order of 5% atomic percent 
[30] a lower value of ^as^^- ~ 8-3 A'^ is obtained. 

For the former case of LaOo.gFo.iFeAs we estimated these corrections due the 
presence of other ions, too. Using the following values for their polarizabilities |52j : 
1.13 A3 for La3+, 2.68 for O^", 3 A^ for F" [50], and again -1 A^ for Fe'^+ [50]. As a 
result we finally arrive at an empirical value of d^s^- = 10.49 A'^. But using again the 
empirical polarizabilities recommended by Shannon and Fischer [18] with q;q2- =1-988 A^ 
and = 1.295 A^ one arrives at aj^^a- = 8.74 A^. Thus, also for this different set the 
As-polarizabilities of two different materials are large and close to each other [19] . Thus, 
the empirical numbers derived from our reflectivity data are in excellent agreement with 
a As ~ 10 A^ suggested by Sawatzky et al. [H], if the ionic polarizations of all ions are 
taken into account. 

Next, LaOFeP could be used to understand the role of polarization effects 
comparing the £00 with those of the La-1111 based arsenic compounds with sufficiently 
higher superconducting and magnetic transition temperatures, except the parent 
compound itself. Since the corresponding fit parameters except the unscreened plasma 
frequency were not shown in reference [37] we are restricted to a qualitative discussion. 
From the refiectivity data shown in figure 1 of reference [37] one might ascribe the 
local minimum in between 4000 and 5000 cm~^, (i.e. around 0.5 eV) to the plasma 
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Figure 7. Total polarizibilitics from equation ([9]) (left). Schcmatical view of the 
"screened" on-site Coulomb repulsion Ud of Fe 3d states (the so called Hubbard U) 
vs. background dielectric constant (right). For the bare case we adopted Uq = 14.7 eV 
[S7] . The stripe denotes the region of Ud estimated by various electon spectroscopy 
methods 



edge but then one would cirri VG cit ci Icir ger or close Soo for LaOFeP compared with 
the arsenic counterpiece. Therefore we suggest that the actual plasma edge is masked 
by an interband transition or by a P-vacancy induced absorption like from the As- 
vacancies in two of the above considered arsenide systems. In fact, for a sequence of 
polarizarion differences oas — cip = 0, 1, 2, 3 we estimate (predict) from equation (9) 
£^oo,p = 14.4, 10.2, 7.7, 6.1, respectively. Then the plasma edge would be near 0.487 eV, 
0.58 eV, 0.67 eV, 0.75 eV, where eoo,As = 11-5 and = 133. 75A^ have been used. In 
addition we have assumed that the polarizabilities of the other ions La^~, 0^~, and Fe^"*" 
remain unchanged, i.e. they cancel out in the difference <5as — ap. A P-vacancy scenario 
might also explain the suppression of an antiferromagnetically ordered state in favour of 
superconductivity as observed in LaOFeP below 6K. If indeed ^oo^p is sufficiently smaller 
than 11 to 15 found here for two arsenide compounds, we would predict also a somewhat 
larger Ud value for LAOFeP. Both predictions can be checked experimentally. 

To get some insight in the accuracy of equation (Q and to check the polarizability 
of the next important most competing ionic component O we considered the cuprates 
La2Cu04 and Li2Cu02 in the same way using in addition 0.2 for the polarizability 
of Cu2+ and 0.0286 A^ for Li+ summing up all ionic polarizabilities one obtains 
atot = 24.28 A^ and 11.31 A^, respectively, (both with two formula units per unit 
cell). As a result with unit cell volumes oi v = 187.89 A^ and v = 98.42 A^ realistic 
£00- values of 4.54 and 3.74, respectively, are obtained from equation Using the 
set of polarizabilities proposed by Shannon and Fischer [18] aQ2^ = 1.988 A'^, and 
Q^Cu2+ = 1.23 A^, we arrive at a^^^ = 10.524 A^ and = 3.434 in the free ionic 
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approximation [51]. The analysis of the loss function -Im[l/£(u;)] for Li2Cu02 within the 
framework of a five-band Cu3(i-02p extended Hubbard model [60] reveals with e^xi = 3.6 
a very close value. The calculated ionic charges Qcu = 1.797e and Qo = — 1.972e are 
also very close to the idealized ionic picture with values of 2e and -2e, respectively. 

In the cuprates the background dielectric constant ^oo is dominated by the 
polarizability of 0^~, i.e. by the 2p-component similarly as in our case by the As 
4p states. Anyhow, in spite of the surprisingly good description in both cases more 
sophisticated calculations beyond the isotropic cubic approximation used in deriving 
equation (Q are nevertheless highly desirable [611 [62] in view of the importance of 
the polarization problem for the elucidation of the mechanism of superconductivity in 
the iron pnictides. In this context we would like to note that our estimated empirical 
values of aAs < 12A^ are probably not large enough to ensure a bipolaronic formation 
at low temperature, i.e. a Bose condensation-like scenario of superconductivity in the 
iron based pnictides such as proposed by Sawatzky et al. [HI [12] This objection is 
also supported by the observation of a nearly linear temperature dependence of the 
upper critical field near instead of an expected pronounced positive curvature generic 
for bipolarons Bc2{T) oc (Tc — T)^'^ [10], the Pauli-limiting behaviour observed in 
some cases [22l [29] as well as the clear observation of Fermi surface related features 
such as the de Haas van Alphen (dHvA) effect in LaOFeP [63] and last but not least 
the numerous angle resolved photoemission experiments (ARPES) [6^ [65] . Anyhow, 
our conclusions about unstable bipolarons should be taken with some cautions in 
view of the approximate character of equation ([9] and the fact that our estimated 
polarizabilities are determined from optical frequencies around the screened plasma 
frequencies whereas the dielectric polarizibilities entering the original Clausius-Mossotti 
equation are determined usually in the range of lkHz-10 MHz and include both ionic and 
electronic components. The extrapolation of our results to that quasi-static limit even 
below typical phonon-frequencies remains an open and interesting question worth to be 
considered in future. We admit that such ionic effects (or lattice effects in general) might 
support the electronic polarization effects and favour finally the formation of bipolarons, 
too. However, even when bipolarons are stabilized it is unclear what will their mass 
? Too heavy bipolarons is a general serious problem for most bipolaronic scenarios to 
explain rather high-Tc values. Another possibility would be also the so-called boson- 
fermion scenario (or s-channel scenario) where real bipolarons occur not as prepaired 
pairs already above Tc but occur only below Tc where they coexist with Cooper pairs 
(see e.g. reference [66] and references therein). 

In order to illustrate the polarization effect on the Coulomb repulsion on Fe sites we 
consider again the relationship between Ua and ^oo employing equation ([9]). According 
to reference [11] the actual on-site Coulomb repulsion Ud in a multiband Hubbard-model 
(for the sake of simplicity we ignore slight differences between the various involved five 
Fe 3(i-orbitals) is given by the difference 



Ud = Uo- 2E, 




New insight into the physics of iron 



19 



where the polarization energy Ep reflects the electric field produced at all surrounding 
polarizable sites due to the additional charge related to the on-site double occupancy and 
Uo is the bare Coulomb repulsion for which we adopt 14.7 eV [67]. In the point charge 
approximation one has Ep = J^c^ji^"^ / Rfj resulting in 2ae'^/Rp^j^g = 0.868q!/[A'^] eV in 
the nearest neighbour (NN) approximation ignoring the infiuence of other ions, where 
RFe,As ~ 2.4 A denotes the NN arsenic-iron distance. Instead of taking into account 
also the dipole-dipole interaction and summing up all contributions within a single 
FeAs layer as it was done in reference [H] in the point charge approximation, we simply 
assume that the polarization energy Ep scales linearly with a where the constant of 
proportionality ui has been determined from the condition to reproduce approximately 
our DMFT based result of ^oo = 5.4 for U = 4 eV given in the previous subsection 4.1. 
Using equations (9) and (11) the resulting phenomenological curve shown in figure 7 
(right) reads 

3v (gpo - 1) 

where v denotes the volume per iron site. This way we arrive at Ud ~ 1.5 eV for 
LaOo.gFo.iFeAs and U ^ 1 eV for BaFe2As2 using the same DMFT based coefficient ui 
also for the latter compound. Thus, our estimated Ud values are in both cases smaller 
than the total band width W ~2 eV. Taking into account the quasi-degeneracy of 
all five Fe 3d orbitals involved one comes to the conclusion that repulsive correlation 
effects should play a minor role except maybe for the magnetic exchange interaction 
affected by the superexchange and the Hund'rule couplings as well. The same reduced 
size of the on-site Coulomb repulsion has been estimated based on various electron 
spectroscopy measurements [68l [69]. Note that residual Coulomb interactions of 
Ueff = 320 to 380 meV and a Hund's rule coupling Jg// = 70 meV have been used 
for instance by Korshunov and Eremin [701 [7T] within a four-band theory for the 
description of thermodynamic properties related to the nesting derived spin density 
wave (antiferromagnetic) ground state and its magnetic excitations of the undoped 
parent compound LAOFeAs and other the iron pnictides (We remind the reader that 
the extended Hubbard model including besides all five Fe 3d also the As Ap orbitals is a 
16-band model.) Without doubt our estimated strongly reduced Ud ~1 eV-values might 
be helpful to understand microscopically the low value of the former, especially, if a 
reduced, but yet sizable, value of the NN intersite Fe-Fe Coulomb interaction V is taken 
into account: Ueff ^ U — V, provided a repulsive \^ > is still obeyed in contrast to an 
expected attractive V < adopted within the bipolaronic scenario of Sawatzky et al, 
where the polarization oas > 12A^ is so strong to change its sign, i.e. V would be over- 
screened. In addition, Ueff is also reduced in the process of mapping down the 16-band 
model "integrating out" the As-4p states [56] exhibiting a moderate hybridization with 
the Fe 3d states. Anyhow, the extremely small value of Jeff compared with expected 0.7 
to 1 eV remains as a puzzle to be resolved. (Since 2 J is given as the difference of the both 
onsite (intra-atomic) U and the intra-atomic-interorbital U' and to first approximation 
both are expected to be reduced in the same way by the polarization of the surrounding 
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ions.). Possibly, the explicit account of an additional direct ferromagnetic exchange 
interaction Kp^ might be helpful to resolve this puzzle. Kp^ ~ 200 meV was introduced 
for the two dimensional corner-shared cuprates by Hybertsen et al. [72] . For instance, in 
almost all edge-shared cuprates such as the above mentioned Li2Cu02 or in Li2ZrCu04 
[71] with Cu-O-Cu bond angles near 90° (just like as the Fe-As-Fe bond angle in the 
iron pnictides considered here) Kp^ amounts about 50 to 80 meV [60], [71]. There Kp^ 
is decisive for the observed ferromagneic NN-exchange interaction being much more 
importantant than the ferromagnetic coupling provided by the Hund's rule coupling Jp 
on the intermediate oxygen O 2p orbitals. The reason for that is that Kp^ affects the 
total exchange already in the second order of the perturbational theory whereas the 
Hunds' rule Jp coupling appears only in fourth order (four small pd transfer intregrals) 

Also a negative [/-scenario as proposed in reference [75] to solve the puzzle of 
the too large calculated magnetic moment in the framework of LDA-|-[/ seems to be 
unlikely. According to our estimate ^oo > 25 would be required (see figure [Tj). Anyhow, 
a detailed comparison with other metalic compounds which exhibit large e-values would 
be interesting also in view of searching for new superconducting compounds and for 
systems with other exotic ground states. 

First estimates of the corresponding polaronic effect result in mass renormalization 
by a factor of 2 in accord with our findings for VLp and those from ARPES data for 
122 superconductors [M] and de Haas van Alphen studies for LaOFeP [63]. Thus, the 
mentioned above strongly reduced Coulomb repulsion Ud can be understood based on 
empirical considerations for the ionic polarizabilities and is in accord with empirical 
findings analyzing various electron spectroscopy data [Ml [Ml [56] . 

5. The el-boson coupling strength and constraints for pairing mechanisms 

Naturally, the most interesting problem of the pnictide superconductivity is the 
elucidation of its mechanism. At the present relatively poor knowledge and still not 
very high sample quality (with respect to impurities and other defects) reflected for 
instance by the high residual resistivities po and an unknown mesoscopic structure due 
to a pronounced tendency to an electronically driven phase separation [76l [77] at least 
for the (122)-family, a comprehensive theoretical description is impossible. In such a 
situation one may adopt various simplified scenarios (with increasing complexity) and 
look for qualitative checks like bipolaron condensation vs. Cooper pair formation, the 
symmetry of the order parameter (s), the coupling strength and the effect of disorder in 
order to discriminate at least a part of the huge amount of already proposed scenarios. 

5.1. One-band Eliashberg analysis 

We start with the simplest possible version of an effective single-band theory where the 
formation of a condensate of Cooper pairs is assumed. (The case of prepaired bipolarons 
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can be excluded by the nearly linear T-dependence of the upper critical field near 
(see Fuchs et al. [22l [29].) Here the mass of the paired charge carriers is affected by 
the virtual retarded interaction with collective excitations (bosons) being present in any 
solid. The corresponding el-el interaction may be attractive or repulsive, supporting or 
weakening a specific type (symmetry) of pairing in addition to the unretarded Coulomb 
interaction which may be also repulsive or attractive under special circumstances |16j . 
Thus, from the standard Eliashberg-theory applied to the case of type-II superconductor 
with weak or intermediately strong coupling (A < 2) [78] one obtains for the penetration 
depth Ai(0) at zero temperature using convenient units 

,Ai(0)[nm] 



v^^p[eV]^^;^^ = v/Wn.(0))(l + A,,,(0))(l + 5), (13) 

where A2,(0)[nm] is the experimental (a, 6)-plane penetration depth extrapolated to 
T = 0. VLp = ^/e^u!p denotes the empirical unscreened plasma frequency shown in 
table 2, ns(0) is the density of electrons in the condensate at T = and n denotes the 
total electron density of all conduction bands which contribute to the unscreened ilp, 
and 6 is the disorder parameter 

e ^ p, y limp _ p. y lirnp ^ limp /. ,n 

2A(0) ~ RT, ^ 5T, ' ^ ' 

which is small but finite even in the quasi-clean limit. R denotes the gap-to Tc ratio 

R = 2A(0)/Tc which in moderately strong coupled superconductors slightly exceeds the 

BCS-value Rbcs = 3.52. Notice that all three factors, the strong coupling correction 

described by the mass enhancement factor Z = 1 -|- A, the disorder D = 1 + 6, as well 

as the condensate occupation parameter = n/ug under the root obey the inequality 

> 1. Let us note that the possibility to extract the e/-6oson-coupling constant from 

a comparison of the high frequency optical response with the renormalized low-energy 

plasma frequency which determines the superfiuid density has been proposed (in the 

clean-limit D = 1 and for the BCS case = 1) already in 1987 by Gurvitch and 

Fiory [79] and again in 1992 by Hirsch and Marsiglio [80]. But to the best of our 

knowledge there were no systematic applications to real superconductors. A first step 

in this direction is presented in subsections 5.1.1 and 5.1.2. Note that equation (fT3ll 

can be generalized taking into account two unusual situations which, however, may be 

present for certain pnictides as discussed below in subsection 5.1.2 in the context of 

LaOFeP [37] in more detail: (i) some residual mass renormalization is still present at 

room temperature in the frequency region of the screened plasma energy ~ 0.4 to 0.5 eV 

i.e. A(u; = Up, T = 300 K) is not very small. This may happen when the spectral density 

of the coupled boson extents to very high-frequencies and (ii) the superconductivity is 

somewhat suppressed due to fiuctuations of a competing neigbouring (e.g. magnetically 

ordered) phase. To take these effects into account we may write: 

1 + Atoi(0) Z-1 Xtot{uJp,T) 

or Xtot{0) = , where r = — - — — — , (15) 



l + \ot{ujp,Ty ' l-r' Atoi(O) 

and 

D = (l + 5)(1 + /), (16) 
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where / > is a phenomenological parameter which measures the fluctuational 
influence of all competing phases. The latter effect occurs beyond the standard mean- 
field approximation and a for a treatment of all instabilities (various superconducting 
phases with different symmetry of the order parameter (s, s^,p, (i- waves), SDW's, 
ferromagnetism, and CDW's, etc.) on equal footing. Since such a corresponding 
comprehensive microscopic theory is not avalaible at present, / has to be treated as 
a phenomenological parameter. 

In classical BCS-type superconductors at T = all quasiparticles participate in 
the superfluid condensate, i.e. = 1. In complex multicomponent superconductors 
like underdoped cuprates the possibihty of < n has been discussed by several 
authors (see e.g. Riifenacht et al. [81] ). In the present case it might also happen 
in the case of coexisting CDW or SDW orderung with superconductivity as observed 
for instance in transition metal dichalcogenides such as NbSe2 or in clean magnetic 
borocarbides (HoNi2B2C). The well-known classical superconductor Nb which much 
stimulated the development of the strong coupling theory provides an almost ideal 
possibility to illustrate the power of equation f|T3|) . Adopting = 1 we may examine 
the case of very clean samples with resistivity ratios RRR up to ~ 2-10^ (0[S1] for which 
a small disorder parameter 6 is expected. In fact, with Al(0) = 31.5 nm [82], 32 nm 
[83] . 34.3 nm [H], 35 nm [85] and the calculated IDA- value of Vtp =9.24 eV [86] or the 
empirical value of 9.2 eV [85] and A = 1.12, = 9.22 K [87] one arrives at 5 = 0.026, 

0. 059 and 0.217, respectively, which point to a very clean-limit regime. In fact, using 
the experimental gap to Tc ratios of 4.1 [HSllHH] and 3.79 [HI] these (5-values correspond 
very small scattering rates of 1.4 K, 3.2 K, and 10.8 K, respectively, only. Based on 
these results we will denote a regime with 5 < 1 as a quasi-clean limit. 

Turning back to LaOo.gFo.iFeAs we adopt also the BCS relation n=ns and 5 ~0.93 
being still in the quasi-clean limit at T = 0. Then, using = 12 [28] and 
Al(0) = 254 nm [5] or 242 nm [89] one estimates Xtot ~ 0.6 and 0.45, respectively, 

1. e. the superconductivity is in a weak coupling regime for our samples with a bulk 
Tc = 26 K (according to the muSR-data [5] while the resistivity yields a slightly larger 
value of 27-28 K). 

The effect of the dielectric background constant renormalized by both the disorder 
and the actual condensate density e* = Soo/DN on the coupling strength is shown in 
figure [HI Note that a substantial impurity scattering absent in the quasi-clean limit 
would further reduce Xtot.. We note once more that our empirical VLp 1.36 eV differs 
from the LDA-prediction of 2.1 eV for the nonmagnetic ground state within the LSDA. 

5.1.1. Comparison with other pnictide superconductors In the present context it is 
interesting to apply our analysis also to other iron based pnictide superconductors for 
which the penetration depth is known. Khasanov et al. [90] reported /zSR measurements 
for the related FeSeo.ss system with a relative low Tc = 8.26 K (probably related to the 
smaller polarizability of Se compared with As). They found from their /xSR-data a 
penetration depth Ai(0) = 406 nm. Adopting a typical plasma frequency of about 
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1.3 to 1.5 eV (compare table 2), one arrives at ZDN ^ 7.16 to 9.53. For a weak 
to intermediately strong el-boson coupling strength Xtot = 0.5 to 1 one is left with a 
relatively large factor DN > 3.58 to 4.76 i.e. either the system under consideration is 
strongly disordered or affected also by fluctuations of a competing phase (see equation 
( fT6l) ). Alternatively this factor is enhanced due-to > 1, i.e. the condensate density 
Us is much smaller than the total electron density n. Finally, Li et al. [31] studying 
Bao.6Ko.4Fe2As2 {T^ = 37 K) reported Ai(0) = 200±8 nm adopting a dirty limit scenario 
in order to describe the observation of a gap in their reflectivity measurements. Then, in 
fact, adopting a relatively large value of fip ^ 1.8 to 2.1 eV such as suggested in reference 
[3T] for the closely related system Bao.55Ko.45Fe2As2 one estimates from equation (fT3l) 
ZDN= 3.28 to 4.47 and there is room for some disorder for weak to intermediately 
strong coupling Z ~ 1.6 to 2 provided A^ = 1. 

The most interesting from the disorder and high-field properties point of view is 
the closely related system given by the arsenic-deficient compound LaOo.gFo.iFeAsi.^ 
with improved superconducting properties near Tc ~ 28 K [22l [29] . The analysis of the 
Pauli limiting (PL) behavior for the closely related LaOo.gFo.iFeAsi.^ system yields a 
similar value A ~ 0.6 to 0.7 as derived from the Eliahberg-theory corrections to the BCS 
(called often "strong-coupling" corrections even for small or moderate values of A ) for 
its PL field -Bp(O) = 102 T [22], [29]. Quite interestingly, the observed kink in the electron 
dispersion in the ARPES data by Wray et al. [91] for the 122 (Sr,Ba)i_^(K,Na),;Fe2As2 
can be described with an effective A ~ 0.6 of the same order as that found for the La- 
1 Ill-system [21] and adopted above for the corresponding rare earth system (see figure 
9) and discussed in more detail there. 




Figure 8. Empirical relation between the total coupling constant from the 
mass enhancement entering the penetration depth Al(0) vs polarizability related 
renormalized dielectric background constant e* . The in-plane penetration depth values 
X1^{0) have been taken from references [5l [89]. 
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In spite of relatively large error bars it is noteworthy that an Uemura-plot type 
relation between the superconducting transition temperature and the squared of the 
empirical unscreened plasma energy VLp seems to hold (see figure 9) although the 1111- 
pnictides occur on another curve than the 122-pnictide superconductors [8]. The 1111- 
pnictides are more or less close to the hole-doped cuprates on an almost linear curve 
but with an additional negativ constant. At present it is unclear what does this 
penomenological result imply? Is there a critical carrier concentration necessary for 
the occurence of superconductivity at variance with the Uemura-plot prediction? Or 
does it mean that for the present crystal structure types (1111, 122, and 111) very low 
carrier phases are simply not realized for chemical reasons and the vanishing Tc near 
fip ~ 1 eV is caused by other reasons not directly related to the carrier density? In the 
La-1 Ill-family the smallest Qp occurs just near the phase boundary to the SDW-phase. 

The 122-systems can be described by another linear relation but also with a negative 
constant. Note that the stoichiometric compounds LAOFeP [37] and NaFeAs[17] occur 
far from these two curves for the 1111 and 122 As-based families Further experimental 
and theoretical work is required to settle these points. A solution might be also 
helpful to understand better the mechanism of superconductivity in iron based pnictide 
superconductors. In the context of a (possibly) somewhat smaller £00 and also a larger 
Ud for LaOFeP, its lower = 6 K as compared with values of 18 K for LiFeAs or 26 K 
for LaOi-ajF^^FeAs might be a relevant hint for the elucidation of the still unknown 
mechanism of superconductivity in Fe pnictides, too. 

The measured in-plane penetration depths of two LiFeAs powder samples with 
superconducting transition temperatures Tc = 16 and 12 K, respectively, amount 
A)L(O) = 195 nm and 244 nm [7j. Since to the best of our knowledge optical data 
are still not available we estimate Qp =1.77 eV using equation (7) and the reported 
lattice constant a= 3.77 A. The isostructural NaFeAs exhibits a larger lattice constant 
a=3.947 and a lower Tc=9 K [17] . 

LAOFeP deserves special attention due to its unusual properties and the relatively 
large amount of experimental single crystal data available for a phenomenological 
comprehensive analysis [33, |92|, [93]. The penetration depth exhibts a quasi-linear 
temperature dependence pointing to unconventional superconductivity with nodes in 
the order parameter [92] which points possibly to a different superconducting phase 
compared with the nodeless FeAs based superconductors. From the measured in-plane 
penetration depth Al(0) = 240 to 250 nm [92] and Qp = 1.85 eV we estimate a large 
value for ZDN ^ 5.06. Ignoring the possibility of > 1 and adopting A^t < 1 as 
suggested by the dHvA-data [63] and ARPES data [95], one is left with D ~ 2.5 to 3.4. 
Since the measurements were performed on very clean single crystals [92] a large disorder 
parameter alone seems to be not very likely to explain the data. Instead we assume a 
sizable contribution of fluctuations measured by our phenomenolgical parameter / (see 
equation f|T6|) ) which is thought to be caused by the presence of at least one competing 
superconducting and/or magnetic phase (s). In this context we note that the possibility 
of enhanced penetration depths due to the vicinity of fluctuating magnetic phases have 
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been assumed very recently also in reference [HI]. Finally, we note that the unsaturated 
temperature dependence of the optical mass renormalization shown in the inset of figure 
3 (panel b) in reference [37] points to a residual \{uj = Up, T = 300K) introduced 
in equation (fT5|) which might be caused by a very broad spectral density (Ehashberg 
function) a{uj)F{uj) with non-negligible high energy contributions. Adopting r(300K) = 
0.25 to 0.33 would be helpful to reduce the differences of mass renormalization as seen in 
optical data (A ~ 0.5) compared with dHvA [63] and ARPES data [95] pointing to A ~ 1. 
The account of r in Z would further increase D and the phenomenological fluctuation 
parameter /. If such an effect would be present also in other pnictide superconductors, 
our estimated coupling constants would be enhanced to Xtot ~ 1- 

Anyhow, whether the quantity Qp is related in fact to the superfiuid density Ug 
remains unclear. Due to some similarity with the original Uemura-plot one might 
conclude that all free charge carriers are in the condensate and the renormalization 
by the el-boson interaction should be rather similar. Then the only free parameter 
within a single band approach is the disorder parameter 6. In this approach the rare 
earth 1111 compounds with the highest Tc- values ~ 40 to 56 K exhibit also the weakest 
disorder. For instance, adopting the same Xtot = 0.61 as for LaOo.gFo.iFeAs discussed 
above we take into account the increase of Qp due to the shortening of a as modelled 
by equation ([7]) and arrive at 6 ~ 0.5. This is a clear indication for the presence of 
multiband effects. Finally, disorder affects mainly that band with the slow particles to 
which the penetration depth A^ is not very sensitive (see below). Finally, we note that 
in the related nickel arsenide LaOo.gFo.iNiAs a coupling constant of A ~ 0.93 has been 
found from a single-band Eliashberg analysis of its thermodynamic properties |35j . 

5.1.2. Comparison with other exotic superconductors Finally, it is interesting to 
compare our results for the Fermi surface averaged total coupling constant Xtot of 
the Fe-pnictides given above with a similar analysis for other more or less exotic 
superconductors where the symmetry of the order parameter may be unconventional 
and in addition to phonons also other bosons such as crystal field excitations and/or 
spin fluctuations may play a crucial role. 

Let us start with extended s-wave superconductors with intermediately strong 
el-ph coupling X^i-ph ~ 1 for which the necessary experimental data is available (a 
corresponding comparison with LDA-based bandstructure calculation will be given 
elsewhere). We start with 

(i) the nonmagnetic borocarbide YNi2B2C (T^ = 15.6 K). According to Widder et 
al. [96] one has VLp = 4.25 eV for the plasma frequency already renormalized by the 
el-ph interaction and Xopt = 1-2 in accord with the Fermi surface averaged A = 1.05 from 
specific heat data [97] . Using the penetration depth Al(0) ~ Xl(T = 3K) = 58 nm from 
the yuSR data by Ohishi et al. [98], one arrives at (1 + 6) = 1.56. Hence, a quasi-clean 
limit value such as 5 ~ 0.7 estimated from the low-temperature optical data of Bommeli 
et al. [99] yields Xtot ~ 0.92 rather close to the value of 1.05 derived from specific heat 
data mentioned above. Similarly using Qp = 6.01 eV and 6 = 0.7 [99] one arrives at 
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Figure 9. Ucmura-plot type relation between tiie superconducting transition 
temperature Tc and the empirical squared unscreened plasma frequency not 
renormalized by the el-boson interaction as derived from optical and ^SR measurements 
for selected Fe-pnictide superconductors. The shown materials are LaOi_2;F3;FeAs 
(1- a;=0.06 [8], 2-x = 0.075 [5] 3- x^O.l [21]) 4- LaOo.gFo.iFeAso.g (present work), 
5- Bao.55Ko.45Fe2As2 [31], 6 - Bao.6Ko.4Fe2As2 [S], 7 - CeOo.84Fi^o.i6FeAs [4], 8- 
NdOo.88Fo.i2FeAs [4] 9 - Smo.82Fo.i8FeAs [89], 10- NdOo.ssFeAs [6] 11 - SmOo.ssFeAs 
[6]. The points 7-11 and LiFeAs are predictions extrapolated from equation ([T]) using 
the actual lattice constant a. 



X = 0.996 even still closer to the specific heat value. Next, let us consider 

(ii) a typical Al5-superconductor, namely, VsSi (Tc = 17 K). Here Al(0) ~ 
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\l{T = 3.8K) = 108 nm derived from /xSR measurements by Sonier et al. |1UU] . The 
experimental plasma frequency of Qpi =2.8 ± 0.2 eV P2| llOlj has been obtained from 
room-temperature infrared data in between 0.16 and 0.62 eV using a pure Drude-model 
relation according to it the slope 1/(1 — ei) vs cu^ is given by {Qp/h)""^. But this 
means that no ^oo has been taken into account. Hence an Qp obtained this way actually 
represents a screened plasma frequency and the real unscreened plasma frequency must 
be larger. For this reason one obtains too small disorder parameters 6 if one one 
adopts for total coupling el-boson constant Xtot = 0.9 obtained from an Eliashberg- 
theory analysis of its thermodynamic properties [102] based on tunneling data for the 
spectral density a'^F{uj) [103] . This procedure is the most accurate way to determine 
the relevant coupling strength for a standard superconductor. Then using the above 
mentioned penetration depth one would have (1 + A)(l + 6) =2.026, 2.35 to 2.7, i.e. 
S = 0.026, 0.24 to 0.42 for the lower bound, the mean value and the upper bound, 
respectively. Thus it is clear why the lower bound is unrealistic. We may refine this way 
Qp = 2.9 ± 0.1. Then one is left with a ratio of 1.15 between Qp^hDA =3.35 eV and the 
refined experimental value. This ratio is a bit smaller than ~ 2 found for the pnictide 
superconductors considered above. Another phonon mediated superconductor with a 
sizable el-ph coupling is 

(iii) MgCNis for which A ~ 1.8 has been found and flp^LDA = 3.17 eV has 
been calculated applying the mentioned above FPLO-code [78]. According to table 
2, and according to our general experience for empirical plasma frequencies, one has 
flp ^ {2/3)Qp^LDA- Thus we adopt flp = 2.1 eV. Using the measured penetration depth 
of At(0) = 231.5 nm |104j one arrives at DZ = 6.07 substituting our empirical strong 
coupling value A = 1.8 we are left with a relative small disorder parameter 5 = 1.17 ~ 1 
in accord with our previous assignment of a quasi-clean limit regime in this compound 

(iv) The Bai_a;Ka;Bi03 (Tc ~ 30 at optimal doping) exhibits a penetration depth 
of A/,(0) = 198.5 nm |105j and A = 1.4 has been estimated by Zhao [105j . With an 
unscreened plasma frequency of Qpi =2.9 to 3 eV |106l 1107] one arrives at relatively 
large values 6 ~ 2.55 to 3.25, where in the last case A ~ Xtr ~1 suggest by Nagata et al. 
from the high-T resistivity data using equation (2) has been suggested. In other words 
Bai_3;K3;Bi03 is a dirty strongly coupled superconductor. 

Turning (v) to the heavy fermion superconductor PrOs4Sni2 (Tc = 1.85 K) and 
to its non-4/ counter part LaOs4Sni2 (Tc = 0.74 K). Their penetration depths amount 
Al(0) = 353.4 nm and 470 nm, respectively. To the best of our knowledge there are 
no optical data available. But theoretical calculation predict very strongly coupled 
superconductivity with A ~ 3.3 in the former case and weakly coupled superconductivity 
with A ~ 0.3 in the last case [108] . Here the strong coupling contribution is expected 
due to nonphononic quadropolar crystal field excitations. If this picture is correct, one 
might predict lower bounds for the unscreened plasma frequency: flp > 1.16 eV in the 
former case and 0.479 eV in the latter case. Finally, turning to cuprates, we consider 
first 
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(vi) the hole (under) doped HTSC YBa2Cu3 06.67 with Ai(0) = 178.4 nm (see figure 
6 of reference [109] ) and f2p ^ 2 eV one arrives at (l + A)(l+5) =3.27, i.e. strong couphng 
Xtot ~ 2 in the quasi-clean limit regime which is in nice accord with the ARPES data of 
the underdoped YBa2Cu306.6 system where A^t = 1-59 has been derived from a kink in 
the dispersion-laws of the quasiparticles |110j and it was described by the authors mainly 
to interaction with the magnetic resonance mode. We close this subsection considering 

(vii) the electron-doped Ndi.85Ceo.i5Cu04_5 where the screened plasma frequency 
derived from the peak of the loss function amounts about 1.06 eV and the plasma edge 
(minimum) in the refiectivity data occurs near 1.3 eV [lllj . The volume of the unit cell 
ammounts v = 188.7 A3 [112] . Using the ionic polarizabilities like those mentioned in 
subsection 4.2 we arrive at an unscreened value of Qp in between 2.2 and 2.7 eV. Using 
equation ([9]) we estimate £00 ~ 4.51 and this way 2.25 eV < ^Ip < 2.76 eV. Anlage et 
al. |H3] measured Al(0) = 105 ± 20A. Whereas Luke et al. found in /iSR measurements 
Al(0) = 230 nm. Then we finally estimate XtotiS^) < 1 for this cuprate with a Tc ~ 20 K, 
only. Our estimated coupling constant Xtot is in accord with a recent theoretical estimate 
of 0.7 by Cappelutti et al. |113] based on the Holstein-t-J model and of 0.8 by Park 
et al. |114j based on an analysis of their ARPES data. Hence, electron-doped cuprates 
exhibit a similar behaviour as LiFeAs and NaFeAs with respect to the relation of and 
the unscreened plasma frequency (see figure 9) and with respect to the el-boson coupling 
constant with La-1111 iron based superconductors. 



5.2. Multiband effects 

Let us consider the simplest case of two effective bands, where for instance one 
band could stand for the two electron bands and the second band for the two hole 
bands predicted by the DFT-LDA calculations mentioned above. Rewriting the total 
unscreened plasma frequency of equation ([6]) as 

nl^t = ^l + ^l^^l{'^+f3'), (17) 

where (3 = 1^2/^1 and the index 'p' has been ommitted, one obtains the following 
straightforward generalization of equation ( |T3i) (compare also the two-band expressions 
for the penetration depth as given in references [1151 1116] ): 
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ZiDiNi 
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where denotes the mass enhancement factor of the band a = 1,2. = 1 -|- Aq, 
with the total coupling constant (i.e. summed up over all interacting el-boson channels) 
including also the interband scattering, Na = na/ua^s is the inverse relative condensate 
occupation number of the band a and Da = 1 + 5q, is the corresponding disorder 
parameter generalizing equation (fT4l) . For coinciding coupling constants and disorder 
parameters, i.e. Ai = A2 and Si = 62, the one-band picture with its weak or moderate 
coupling strength discussed above is reproduced (see figures 8 and [10]). Upon weakening 
(strengthening) the coupling in one band say, in band 1, stronger (weaker) couplings are 
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Figure 10. Constraint for the el-boson coupling constants Ai and A2 for various 
individual Drude plasma frequencies measured by the parameter /3 = ^2/^1 using the 
experimental value of the penetration depth Xl{0) =254 nm [5] and the empirical total 
unscreened plasma frequency flp =1.37 eV [Hj for LaOo.gFo.iFcAs (left) the same for 
enhanced disorder in the "slow" band 2 (right). 



introduced in band 2. Naturally, the band with the larger coupling exhibits the larger 
gap, if the corresponding couplings support the actual pairing state which, however, has 
not been specified in our approach. From that constraint there is, in principle, taking 
into account the present poor knowledge of parameters no practical bound for A2, i.e. for 
the band with the slower quasi-particles, whereas the coupling constant Ai for the band 
with the faster charge carriers is restricted from below and from above as well. At the 
upper bound A2 = 0, whereas it formally diverges at the lower bound, provided A™" > 0. 
In the special case when <^ 1, i.e. of a light and a heavy band, the fast electrons may 
mask the heavy electrons. For the iron pnictides under consideration the hole bands 
do exhibit somewhat smaller Fermi velocities (and also plasma energies) according to 
band structure calculations. Therefore one might ascribe our effective band 2 to these 
hole pockets. In general, we regard equation ( fT8|) as potentially helpful to identify 
the location of the larger and the smaller gap values in addition to other methods. A 
strongly coupled superconducting heavy subgroup which would be responsible for very 
high upper critical fields might be this way hidden by weakly coupled fast electrons. 



6. Conclusion 



We performed a combined theoretical and experimental study of the screened and 
unscreened in-plane plasma frequencies of iron based pnictide superconductors. The 
obtained moderate renormalization of the corresponding optical masses is similar the 
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one found frequently for other transition metals and is probably a normal Fermi 
liquid effect which gives strong support for the relevance of the electronic structure 
as calculated by various density functional based methods. Evidence for the presence 
of weak correlation effects can be deduced also from the obtained large background 
dielectric constants, exceeding 11, which is ascribed to the large polarizability uas- A 
comparison with zero temperature penetration depth data from muon spin rotation data 
provides a new constraint for the total coupling strength of the el-boson interaction, the 
strength of impurity scattering rates and the density of quasi-particles involved in the 
superconducting condensate. This information yields also reasonable estimates for other 
more or less exotic superconductors. 
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